Local origins of volume fraction fluctuations in dense granular materials 
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Fluctuations of the local volume fraction within granular materials have previously been observed to decrease 
as the system approaches jamming. We experimentally examine the role of boundary conditions and inter- 
particle friction /x on this relationship for a dense granular material of bidisperse particles driven under either 
constant volume or constant pressure. Using a radical Voronoi" tessellation, we find the variance of the local 
volume fraction monotonically decreases as the system becomes more dense, independent of boundary condi- 
tion and /i. We examine the universality and origins of this trend using experiments and the recent granocentric 
model fl 2 1, modified to draw particle locations from an arbitrary distribution V(s) of neighbor distances s. 
The mean and variance of the observed V(s) are described by a single length scale controlled by 0. Through 
the granocentric model, we observe that diverse functional forms of V{s) all produce the trend of decreasing 
fluctuations, but only the experimentally-observed V(s) provides quantitative agreement with the measured 
fluctuations. Thus, we find that both V{s) and ^(0) encode similar information about the ensemble of observed 
packings, and are connected to each other by the local granocentric model. 



I. INTRODUCTION 

Recent measurements of granular systems show that 
densely packed aggregates exhibit smaller fluctuations in their 
local volume fraction </) than more loosely-packed systems. 
This has been observed for several static systems (3H21, as 
well as for a dense driven granular system O where particles 
with different frictional properties each exhibited the same 
quantitative relationship. Measurements of fluctuations in the 
global volume fraction $ (or total volume) also demonstrate a 
similar trend in the fluctuations around a steady state value 1 7~ 
[TOl . However, in several cases, the shape of the relationship 
between ^ and its variance did not monotonically decrease 
(SKTO I. This bulk behavior may be related to the onset of co- 
operative effects and a possible phase transition ifTTIl . In con- 
trast, for local measurements of (j) of static 1 3 - 5 1 and dynamic 
161 packings, the decrease in the variance was monotonic in 
(p. The observations of this trend span experiment and sim- 
ulation, various preparation protocols and particle properties, 
and both two and three dimensions, suggesting that a universal 
explanation might underly the observation. 

For dense systems, the decrease in the fluctuations of global 
$ (or local (j)) is suggestive of a decreasing number of valid 
configurations as the system approaches random close pack- 
ing. Heuristically, this can be understood by considering six 
nearest neighbor particles arranged in a ring surrounding a 
central particle. As the size of this ring shrinks (corresponding 
to increasing (j) locally), the number of possible configurations 
for the neighbors decreases until there is a only single config- 
uration possible in a hexagonally close packed state. 

One framework in which to describe the global fluctuations 
in the volume fraction is the Edwards approach to the sta- 
tistical mechanics of static granular systems |12|. The en- 
tropy S{V) = k\ogfl{V) increases with the number of me- 
chanically valid, static configurations ft, which is a function 
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of the system's volume V for a constant number of parti- 
cles. The change in this entropy as a function of V provides 
a temperature-like quantity, compactivity X = dV/dS, for 
which X Oas<l> ^rcp (random close packing) and 
X ^ oo as ^ ^ ^RLP (random loose packing). The vari- 
ance in either F or <l> is associated with the compactivity- 
analog of specific heat. It remains an open question how to 
connect local, statistical measurements of to a 'thermody- 
namics' of the bulk system for jammed systems |[l3l[T4l. It 
is even less clear how one might apply such descriptions for 
dynamic or even slightly unjammed configurations, as in the 
experiments presented here. 

On the particle-scale, the local volume fraction (j) is defined 
as the ratio of the volume occupied by the particle to the total 
locally-available space. One method for partitioning space is 
the radical Voronoi tessellation (also known as Laguerre cells 
or power diagrams), which constructs cells according to the 
locations and radii of the neighboring particles |[T5]|T6l. Each 
Voronoi cell contains a single particle i, with the boundaries 
of the cell enclosing the set of points for which the distance di 
to the particle i satisfies di < -\- — r\ (k are the indices 
of all other particles in the packing). This tessellation tiles all 
space, with one cell for each particle; neighboring particles 
are defined as those having cells which share an edge. The use 
of the radical tessellation is desirable for dense polydisperse 
systems in order to ensure (j) < 1. This choice is not unique, 
and alternate methods for partitioning space at the particle- 
scale have also been utilized to similar effect OJIMI ,15, 17 1. 
Given a complete set of cells which tile the volume, we can 
define a local volume fraction (j)i = Vi/Vi, where Vi = nrf 
is the volume of the i^^ particle and Vi is the volume of its 
Voronoi cell. The mean and variance of this distribution are 
denoted (j) and = ((^ — ^)^), respectively. A bar over a 
local variable indicates a mean over all particles. 

For spheres in two or three dimensions, hexagonal close 
packed order provides the densest packing with a volume 
fraction of = i\ j \fVl ^ 0.91 in two dimensions and 
(\)\^ = 7r/Vl8 ^ 0.74 in three dimensions (TSl . For both 
of these ordered packings, only a single local configuration 
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is possible = cj) = const.) and the variance of the dis- 
tribution V{(j)) is thus cr^ = 0. In disordered systems, a 
jamming transition occurs at global volume fraction and 
mean coordination number zj which obey the relationship 
z — zj oc — <l>j)" for packings above the jamming transi- 
tion (T9II20I. In the presence of disorder, static packings ex- 
hibit local variations in both z and (/). As a result, increasing 
the volume fraction and coordination number of the aggregate 
decreases the translational and rotational randomness as well 
as the anisotropy |5, 21 1. Correspondingly, is one measure 
of how much disorder is present in a system, and has the im- 
portant advantage of not requiring knowledge of whether or 
not two neighboring particles are in mechanical contact. This 
makes it an experimentally- tractable state variable. 

Due to the prevalence of the trend of decreasing with 
increasing (j) in multiple jammed and unjammed experiments 
and simulations, we investigate how it arises. We perform 
experiments on an unjammed, driven system to probe the 
robustness of this trend with regard to boundary condition 
(constant pressure and constant volume) and the frictional 
properties of the particles (inter-particle friction coefficients 
/i = 0.04, 0.50, and 0.85). We measure (j) for individual 
particles, and observe that decreases linearly with increas- 
ing (j), independent of boundary condition or /i. As this trend 
is quite similar to observations in jammed systems, we sug- 
gest that geometry plays an important role, rather than driving. 
Thus, we examine the origins of this trend via a generaliza- 
tion of the granocentric model (TJEll, whereby we introduce 
randomness through the nearest neighbor distance distribution 
V{s). We find that while various models ofV{s) all produce 
the trend of decreasing fluctuations, only the experimentally- 
observed V{s) provides quantitative agreement with the mea- 
sured (f) fluctuations. 



II. EXPERIMENT 

The experiments are performed on a single layer of parti- 
cles which are supported by an air hockey table from below 
and driven to rearrange by an array of sixty bumpers which 
form the perimeter (see Fig. [T^, the same apparatus as |6 1). 
The particles are a bidisperse mixture of particles with large 
= 43 mm and small rs = 29 mm radii with masses 
rriL = 8.4 g and ms = 3.4 g respectively. The ratio of 
the number of particles is fixed at 1 large particle to every 
2 small particles. The particles are prepared to have one of 
three inter-particle frictional coefficients: /ii = 0.04 (PTFE 
wrapped), /i2 = 0.50 (bare polystyrene), or /is = 0.85 (rub- 
ber wrapped). The air jets provide nearly frictionless contact 
between the particles and the base, and the table is leveled 
so that particles do not drift to one side in the absence of 
bumper driving. While particles not experiencing collisions 
with either the bumpers or neighboring particles can still drift 
slightly due to the air jets and/or local heterogeneities in the 
table, these velocities are small compared with the dynamics 
induced by the bumpers. 

Each experimental run consists of at least 10^ configura- 
tions captured by a camera mounted above the surface of the 
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FIG. 1: Color online, (a) Photograph (left) and schematic (right) 
of the apparatus, showing confining wall in constant pressure (CP) 
configuration, with weights m suspended from a pulley via a mono- 
filament line. Constant volume (CV) is obtained by fixing the wall 
to the surface of the table, (b) Equation of state (/)(P) for CP exper- 
iments at /ii (+), /i2 (o) and /is (x) on linear axes and inset with 
semilog axes. 



table. Images are collected every 2 to 5 seconds, according 
to the $ -dependent dynamical timescale of the system. The 
configurations are generated by the agitation from the perime- 
ter of the packing by bumpers of width ~ 3.2r5'. A pair of 
bumpers at the same position on opposite walls are simultane- 
ously triggered so that no net torque is exerted on the system. 
Every 0.1 s, two pairs of bumpers are randomly selected and 
fired, maintaining ongoing dynamics in the aggregate. 

The bumpers generate evolving dense particle configura- 
tions at global volume fractions within the range 0.71 < ^ < 
0.81. For reference, we previously measured the static ran- 
dom loose packing for these particles tohe ^rlp ~ 0.81 for 
M2,3 lEl. Experiments and simulations 1 19^, 22", "23 1 with simi- 
lar bidispersity have observed random close packing ^rcp ^ 
0.84. The particle dynamics, driven by the bumpers at the 
perimeter, are caged at short time- scales and diffusive at long 
time-scales |6|. As global ^ approaches jamming, the dy- 
namics slow down sharply, as when the glass transition is ap- 
proached in thermal systems. The system is also well-mixing, 
as measured by the braiding factor 1251 of the trajectories 
growing exponentially in time 1241 . Thus, this apparatus is 
well- suited for generating a large number of sterically- valid 
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(non-overlapping) configurations. 

The boundary condition is determined by a movable wall 
on one side of an approximately 1 m x 1 m region along with 
bumpers on the other three sides of the aggregate. The wall, 
of mass 95 ms, extends the width of the table and can be con- 
figured to provide either constant pressure (CP) or constant 
volume (CV) boundary conditions. For CP conditions, shown 
in Fig.[T] the wall functions like a piston and is pulled towards 
the aggregate by a weight of mass m suspended from a low- 
friction pulley and a mono- filament line. For a fixed number 
of particles N = 186, we perform experiments for a range 
of scaled pressures P = m/rus from P = 0.17 (0.58 g) to 
2.41 (8.14 g). This driven granular aggregate behaves as a 
compressible fluid, with ^(P) shown in Fig.jlj). For CV con- 
ditions, the wall is fixed to the table so that the particles are 
confined within an approximately square region. We vary the 
number of particles N from 180 to 204 (altering the global $) 
while keeping the 2:1 (smallilarge) ratio for all TV. 

From each image, we extract the center and radius of each 
particle, and perform our analysis on only the inner 20% 
(about 40 particles), which reduces ordering effects due to the 
boundary |26|. From the particle positions, we calculate (j) us- 
ing the Voro++ |[27ll implementation of the radical Voronoi' 
tessellation. These local measurements allow us to consider 
the probability density function (PDF) V{(j)), as well as its 
mean (j) and variance a^, as a function of the three values 
of fi and two boundary conditions. In addition, we record 
the neighbors for each particle together with the inter-particle 
distance s which separates the edges of the two particles. The 
distribution V{s) will provide a key input to the granocentric 
model. 



III. RESULTS 

As has been previously observed for static granular media, 
the 7^((/>) functions arise from the inverse of a gamma-like 
function of the free volume (161 [28l. Remarkably, for experi- 
ments with similar (j), but different boundary conditions and /i, 
we find the to have similar mean, variance, and shape, as 
shown in Fig. [2p. This observation typically holds at other (j) 
as well. As in 1 16|, the volume distribution is bimodal because 
large particles occupy a larger volume on average. Across all 
experimental runs, we observe a one-to-one relationship be- 
tween the mean and variance for over the entire range of 
global ^, three values of /i, and two different boundary con- 
ditions. Taking the moments ^ and from each distribu- 
tion, we observe a single, linearly-decreasing trend, shown in 
Fig.|2]3. The linear fit gives the intercept = 0.842±0.002 for 
cr^ = 0, which is close to ^rcp- Similar universality in the 
relationship {(j)) was observed in the same apparatus using a 
different technique for measuring </) and |6|. In addition, a 
similar trend, perhaps with discontinuities or changes in slope, 
was previously observed in experimental and numerical three 
dimensional, monodisperse, jammed packings fS^-^. 

The ^-distribution is sensitive to the location of neighbors, 
where each neighbor determines the boundary of one side of 
the Voronoi cell. The probability distribution of s is measured 
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FIG. 2: Color online, (a) Representative ^(0) for three experi- 
mental runs at ~ 0.780: + for CP and /i2, • for CP and /xs, 
and X for CV and /i2. For the run at CP and /xs, dashed and dot- 
ted lines show for large and small particles, respectively. 
All 7^(0), whether or not they distinguish large and small parti- 
cles, are normalized by the total number of measurements, so that 
-^(0) = Vl((I)) + (b) Mean and variance of V((l)) 
measured from individual Voronoi' cells. Each point is for a single 
experimental run with > 10^ configurations. Inter-particle friction 
is denoted by shape (/ii, o; /i2, V; /is, □) and boundary condition 
by open or filled markers (CP, open; CV, filled). The line is a least 
squares fit with ^-intercept at = 0.842 zb 0.002. 



from particle edge to particle edge (see Fig. |3^) and shown 
for two example runs at different (j) in Fig. |3p. In our experi- 
ment, the shape of V^^^{s) is similar in all our data regard- 
less of boundary condition or /i. After a peak in V^^^{s) 
near Smax = 0.05 (1.8 pixels), the probability falls to- 
wards a flat value until a knee at s ^ r^, after which it starts 
to fall off exponentially with a decay set by s. A closer ex- 
amination of V^^^{s) reveals that the location of Smax is af- 
fected by whether the neighboring particles are large or small. 
For large-large pairs, Smax ^ 0.075 (0.051 r^) and for 
small-small pairs Smax ^ 0.026 r^, while the knee remains 
at 5 = rs and the exponential tail of the distribution re- 
mains unchanged. These values are independent of (j). For 
simplicity, we ignore the distinction between large and small 
particles for V^^^{s) with little loss in accuracy of the model 
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FIG. 3: Color online, (a) Schematic of the central particle and one 
neighbor with radii Vc and r, respectively. The shortest distance be- 
tween the edges of the two particles is s. The shaded region is the 
contribution from this neighbor to the total radical Voronoi' vol- 
ume, (b) jy^^P (s) (dotted) measured for two experimental runs at CV 
and /is, with = 0.747 (thick) and = 0.792 (thin). For compari- 
son, exponential distribution V^{s) — ^e~^^^ (dashed) with A = s 
and ^-function distribution [s] — 5{s — ~s) (solid), (c) Scaled 
mean (A) and standard deviation (x) of s, for all BC and /i, as a 
function of the local volume fraction 0. 



results. (This choice may be inappropriate in highly poly dis- 
perse systems.) In spite of the peculiar shape of the observed 
distribution, both the mean and standard deviation of s appear 
to be smoothly set by ^, as shown in Fig. [s]:. This suggests 
that a single length scale controls the distribution. 



IV. MODEL 

Due to the universality of the cr^(^) trend with regard to 
boundary condition and inter-particle friction, as well as the 
observation of a similar relationship present for other dimen- 
sionality, polydispersity, and protocol El |4l |6l [TOl, we seek 
a geometric explanation for the relationship. We choose as 
our starting point the recent granocentric model |1 , 2 | which 
considers the inherently local origins of the volume fraction. 
The model takes the finite amount of angular space available 
in the vicinity of a single particle, and considers a random 
walk which fills this available space with particles drawn from 
a known distribution. We explore suitable extensions to the 
model and examine the implications for the universal trend in 
shown in Fig.|2]3. 

Three parameters govern the model: the size distribution of 
the particles V{r), maximum available angular space 9c (27r 
in two dimensions, 47r in three dimensions), and fraction of n 
neighbors in mechanical contact, Pz = z/n. The z contacts 
are those which provide mechanical stability for the central 
particle, and contribute to the z = 2d condition for isostaticity 
in d dimensions. For measurements of the jammed emulsions 
which formed the inspiration for the model, Pz ~ 0.4 was 
observed and the constant separation s between non-contact 
neighbors the free parameter used to fit the predicted to 
experiment Q. The model invokes randomness through both 
the radius distribution and Pz, and finds 7^((/>) in quantitative 
agreement with experimental measurements. 

To apply the granocentric model to un-jammed systems 
< ^j) such as this one, we consider several modifica- 
tions to the inclusion of V{s). The original model draws s 
from a binomial distribution containing 5 = and a tunable 
s = const, with probability Pz and 1 — Pz, respectively. As 
shown in Fig (Sj), a much wider distribution of s is observed 
in our un-jammed and driven system. Examining V{s) for 
5 0, we observe that mechanical contacts are rare; an upper 
bound of = 0.04 is set by threshold of the image resolution. 
Therefore, we do not treat contacts and neighbors separately. 

Following the original formulation of the granocentric 
model |1 , 2 1, there is a maximum angular space 6c available 
around a central particle with radius Tc which can be occu- 
pied by neighbors. Each neighbor of radius rj with its edge Sj 
away from the edge of the central particle occupies an amount 
of space 

Qj = 2 arcsin (1) 

Tj + Tc + Sj 

which provides a theoretical range < &j < tt; in the experi- 
ments, only the range O.IItt < Qj < O.TGtt is observed. For a 
collection of n randomly- selected neighbors, the total angular 
space occupied is Yli ^j{^ci ^ji^j)^ two-dimensions 
this sum must be less than 9c = 27r. Each neighbor con- 
tributes V* to the Voronoi cell Vi of the central particle, shown 
as the shaded region in Fig. [3^ and given by 
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To compare with the observed V{(j)), we perform a Monte 
Carlo simulation which draws particles from the 2:1 bidis- 
perse size distribution and sj from a specified V{s). For 
each Tc, neighbors are sequentially selected at random from 
these two distributions. For each neighbor, we calculate the 
©-contribution according to Eq. [T] The random process con- 
tinues until Xli^^ > ^TT, at which point insufficient an- 
gular spaces available for the last randomly- selected parti- 
cle. Only the n neighbors are retained and used to calculated 
Vi = J2i the central particle. This process is repeated 

for 10^ different seeds in order to obtain a distribution of local 
volume fractions (pi = nrl/V* . 

In calculating V* , we make one additional adjustment to 
account for the fact that the rejection of the n + 1 neighbor 
leaves behind a neighbor-less gap of size Qex = 2tt — Qj . 
Failure to account for this gap leads to an overestimation of cj), 
which we correct by apportioning Og^ among the neighbors 
in proportion to the angular space they already occupy. The 
adjusted angle occupied by each neighbor becomes 
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so that 9c = Yli ®j = ^TT. Using this adjusted value, Eq.|2] 
becomes 



tan 
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and (/) = 7rr^ / Vj provides a better model of the local 
Voronoi' volume, as the angular space surrounding the parti- 
cle is now completely occupied by neighbors. Note that the 
construction of the boundary at the half-distance (Tc + sj / 2) 
between the edges of the particles does not result in a cell 
with realistic Voronoi shape. This is also true for a boundary 
drawn at a more Voronoi'-like distance ^((^c + + -^j)^ + 
rl — r'j)/{rc -\- Vj -\- Sj) from the central particle. In either 
case, we find that the model reproduces the observed 4>, but the 
half-distance construction quantitatively predicts cr^(^) better 
than the Voronoi-like construction. Therefore, we use the half- 
distance construction for the granocentric model in the results 
that follow. 

We start from four different V{s) distributions: (1) 
the experimentally-measured V^^^{s)\ (2) pexp,cutf^g^ _ 
pexpf^g < rs), the experimentally-measured distribution 
without the low-probability knee; (3) an exponential distri- 
bution V^{s) = je~^^^ with A = s and (4) a delta function 
V^{s) = S{s — s). The effect of the low probability but large- 
s tail of 7^^^^ (5) is illustrated by p^^p^^^t(s). While V^{s) 
and V^{s) may not be physically realistic, these artificial dis- 
tributions are chosen to show the strong dependence of V{(j)) 
on the functional form ofV{s). Each of the four distributions 
is shown in Fig. [3}) for comparison. In order to examine the 
V^{(t)) which arise from these four P{s), a Monte Carlo pro- 
cess draws a random Vc to form the basis for the i^^ cell of 
the distribution. Around that central particle, sequential ran- 
dom neighbors of size Vj are placed at random separations Sj 
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FIG. 4: (a,b,c) Measured ^(0) (dashed line) and calculated gra- 
nocentric prediction {(j)) (o). The granocentric prediction is cal- 
culated using (a,b) 7^^^^(s) from ¥igJ^ and (c) V'^^^^s) for a sim- 
ulated packing with = 10^ particles shown in (d). 



until the available angular space is used up; this collection of 
neighbors provides a value (pi for that cell. This same Monte 
Carlo method is repeated to generate 10^ values of (pi and 
thereby compute V^{(j)) for of the four V{s) for each experi- 
mental run. Two advantages of a Monte Carlo simulation over 
the semi-analytical techniques of |1 , 2 | are that it permits the 
use of the experimentally-measured V^^^{s) as an input to the 
model, in addition to the ability to redistribute Qex among the 
randomly selected neighbors. 



V. COMPARISON 

Using V^^^{s) as the input to the granocentric model, as 
described above, provides a prediction for V^{(p) which is 
in quantitative agreement with the experimental results. This 
comparison is shown in Fig. |4^b, for the same two runs as 
in Fig. [Sj). As expected, the width of the distribution nar- 
rows with increasing (p towards the random close packed limit. 
In Fig[4j all three panels compare the V{(p) observed in the 
experiment and simulation with the granocentric prediction 
V^{(p). In Fig. [4^bc, we observe quantitative agreement in 
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FIG. 5: Color online, (a) Mean and variance of ^(0) for all experi- 
ments (■, same data as Fig.[2|, compared to granocentric predictions 
drawn from several different s distributions: the experimentally- 
measured distribution V^^^{s) (•), the s < rs portion of the 
experimentally-measured distribution -p^^p^^^^ (^g^ (_^)^ exponential 
distribution V^{s) = je~^^^ (A) where A = Sexp, and V\s) = 
S{s — Sexp) (♦)• The filled markers represent distribution in the ex- 
perimentally observed range of 0.154 < Sexp/rs < 0.339. The 
open markers, and A are V^{(j)) and V^[<^) for s < Sexp. (b) The 
mean number of neighbors versus cj) for V^{s) and V^{s), with the 
same markers as in (a). 



the peak and shape of P(0). However, for denser packings 
(0 > 0.75), the model is systematically high on the low-(/) 
side of 7^(0). 

This comparison is performed on a small number of parti- 
cles over many configurations. To test the model for a larger 
number of particles in a single configuration with a different 
preparation protocol, we use the YADE discrete element model 
[29 , 30 1 with 10^ particles, the same bidispersity as the experi- 
ment, and inter-particle friction coefficient fi2 • The simulation 
is similar to the Lubachevsky-Stillinger algorithm [31 1, where 
each particle starts as a point with r = and is grown linearly 
in time proportional to the desired radius with a damping co- 
efficient of 7^ = 0.3. For the purposes of comparison with ex- 
perimental data (which is not in a jammed configuration) we 
end the inflation algorithm when the volume fraction reaches 
(j) = 0.742. Using the particle positions and sizes, we repeat 
the same (j) and s measurements as for the experiments. 

The measured is shown in Fig [4a and exhibits a 

similar shape to the V^^^{s) shown in Fig.^ with an expo- 
nential decay for r > rs- Note that the peak in V^'^^{s) is 
located at s = 0, whereas the peak in V^^^{s) is finite (but 
small). The difference may be due to the differing preparation 
protocols. In the simulations, as the radii of particles grow, 
the overall energy in the system would increase unless kept 
finite by viscous damping. Thus, the dynamics of the system 
are slowed as the desired global <l> is reached. Since Fig. [4]l 
measures the final result of the simulation once all particle- 



overlaps are eliminated, we ensure 5 > 0; due to the slow 
dynamics, most contacts remain at 5 = 0. In contrast, the 
configurations in the experiment arise through collisions and 
the configuration 5 = is less likely. 

Finally, we are able to compare how well the model can ex- 
plain the experimental results shown in Fig. The results 
are shown in Fig. [5^, where the black squares are all of the 
experimentally-measured cr^(^) from in Fig. [2] We observe 
that of the four proposed distributions of s, the full V^^^{s) 
produces the best agreement and is able to capture not only 
the decreasing linear trend, but nearly the correct quantitative 
values. We find that the exponentially -rare neighbors with 
s > rs provide an important contribution to both the mean 
and variance: when they are removed from V^^^{s) to form 
pexp,cutf^g^^ the moments of 7^(0) are less-accurately repro- 
duced. 

For the two heuristic models of 7^(5), we are also able to 
produce the trend of decreasing variance, but without quan- 
titative agreement with experiments. The exponential distri- 
bution V^{s) is a function of a single free parameter X = s, 
where the mean and variance are equal to s and 5^, respec- 
tively. By smoothly varying A over a range of values consis- 
tent with V^^^{s), we obtain systematically lower variance in 
(j), indicating that the particular shape of 7^(5) is important for 
quantitative agreement. Even for the constant s provided by 
V^{s) (likely more consistent with jammed systems than the 
driven ones described here), a linear relationship for {(j)) re- 
mains, although with significantly less variance than observed 
in the experiments. 

Using V^{s) causes cr^(^) to become discontinuous, unlike 
experimental observations. This can be understood as arising 
from the low degree of polydispersity in the system. In FigjSj), 
the discontinuities in cr^(^) and the average number of neigh- 
bors, n, occur at the same values of (j), which is controlled by 
s. As s increases ((/) decreases), the number of neighbors is 
approximately constant until Sex is greater than the mean 6, 
at which point there is on average room for one more neigh- 
bor which appears as a step in both plots. For distributions of 
V{s) with sufficient variance, cr^(0) and n are observed to be 
continuous. 

Although the experiment does not probe > 0.8 due to 
lengthening timescales |6|, we can explore V^{s) and V^{s) 
for larger (j) (smaller s). These points are shown as the open 
markers in Fig.js] In the limit 5 ^ 0, approaches the same 
non-zero value for both models. For comparison, a very loose 
arrangement of particles prepared to have (j) < 0.6 by the same 
Lubachevsky-Stillinger algorithm used in Fig. [4]: exhibits an 
increase in cr^(0) through a maximum near </) ^ 0.5. Com- 
paring the model V^^cj)) with V^'^^^cj)) for low 0, the model 
underestimates the variance even though ^ is calculated with 
good agreement. 

Thus, we find that the shape and moments oiV{s) strongly 
affect the prediction of the local distribution of (j). Very narrow 
distributions such as V^{s), or distributions which resemble 
only a portion of the experimentally measured distribution of 
5, failed to predict our experimentally measured (j). Nonethe- 
less for dense granular systems, all V{s) produced a mono- 
tonic decreasing relationship between and (j), suggesting 
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that this trend is quite robust. 

VI. DISCUSSION 

Independent of boundary condition and inter-particle fric- 
tion /i, we have found that as the mean local volume fraction (j) 
increases, the variance monotonically decreases in our driven 
granular system. This result is reminiscent of similar trends 
observed in static granular systems |l3H5|[7H2l. in spite of the 
special nature of jammed systems. In particular, many proper- 
ties of static granular systems are due to the inter-particle fric- 
tion. Increasing /i in jammed systems decreases the required 
number of contacts from z = 2d for /j. = to z = d -\- 1 
for /i = oo. In a model of local mechanical stability by Sre- 
bro and Levine |32|, increasing /i increases both the tangen- 
tial force and the maximum angle for which two particles can 
be mechanically stable. This would presumably lead to /i- 
dependence in the distribution of local (j). However, in driven 
systems the number of configurations may well be indepen- 
dent of /i as there is no constraint on z. In unjammed systems, 
the lack of /i-dependence on the local ^-distribution also sug- 
gests that dissipation does not play an important role in de- 
termining local (j). However, prior observations in the same 
apparatus 1 6| indicate that driven, equilibrating subsystems in 
fact have ensembles for which Q = which appears to 

be in contrast with this idea. 

The lack of dependence on boundary condition is also sur- 
prising, particularly given two possible sources of anisotropy 
in the system: only three of the four walls provide driving in 
both CP and CV conditions, and the piston introduces a com- 
pressive force in the case of CP. Nonetheless, little anisotropy 
was observed in the angular distribution of neighbors: the 
bond angle order parameter Qe L21J was indistinguishable for 
CP and CV systems at a the same (j). Additionally, we tested 
V{s) for angular dependence with respect to a reference vec- 
tor in the lab frame and found that V{s) was rotationally sym- 
metric for both boundary conditions. 

In global (rather than local) measurements of system vol- 
ume V, Schroter et al. 1 8 1 observed a non-monotonic relation- 
ship for static aggregates prepared through sedimentation of 
a fluidized bed at different flow rates. There, fell from 
0i?LP towards a transition point, then again rose on approach 
to (j)Rc p • This rise was attributed to cooperative effects within 
a finite number of statistically-independent regions, and may 
be related to the presence of spatial correlations during the 
approach to (I)rcp ED- Aste and Di Matteo |3| compared 
the standard deviation as function of (j) for a wide variety of 
experiments and simulations, and found non-monotonic be- 



havior for the aggregated data across the transition from un- 
jammed to jammed configurations, although not for individ- 
ual controlled experiments/simulations. Fig. [5] suggests the 
interpretation that different experiments/simulations produce 
different V{s) and thus fall on different cr^(0) curves, each of 
which has its own monotonic relation. 

VII. CONCLUSION 

We find a geometric explanation for the universally- 
observed trend of decreasing fluctuations with increasing vol- 
ume fraction. By generalizing the granocentric model [ iJEl 
to take the neighbor distribution V{s) as the single input to 
the model, we find that the variance of local (j) measurements 
exhibits a smoothly decreasing trend as long as the first two 
moments of V{s) are also smoothly decreasing. Therefore, 
it is not surprising that cr^((/)) has been observed to decrease 
with (j) for a variety of preparation protocols and particle prop- 
erties in both two and three dimensions. When V{s) is chosen 
to be the experimentally-observed distribution, the granocen- 
tric model provides quantitative agreement with the observed 
shape of V{(j)) and cr^(0) without reference to any spatial cor- 
relations in local (j) due to cooperative effects. Interestingly, 
this suggests that V{s) and P(0) both encode similar infor- 
mation about the distribution of free volume, with the first two 
moments set by (j). While 7^(0) (and the related free-volume 
distribution) has been well-studied in granular systems, V{s) 
has the advantage of being closely related to the radial distri- 
bution function, with the caveat of only considering neighbors 
instead of all particles. 

In conclusion, the decreasing relationship cr^(^) with re- 
spect to increasing reveals the key role played by the con- 
strained availability of angular space in the vicinity of a sin- 
gle particle. This constraint holds regardless of whether the 
system is static or dynamic, jammed or unjammed, mono- or 
polydisperse, or two or three dimensional, and requires no in- 
formation about the force chains or dynamics of the system. 
Therefore, this relationship provides a way of examining the 
distribution of free volume, and perhaps ultimately the full 
ensemble of valid configurations, in both jammed and driven 
granular materials. 
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